Oscillation effects on high-energy neutrino fluxes from astrophysical hidden sources 
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High-energy neutrinos are expected to be produced in a vareity of astrophysical sources as well 
as in optically thick hidden sources. We explore the matter-induced oscillation effects on emitted 
neutrino fluxes of three different flavors from the latter class. We use the ratio of electron and tau 
induced showers to muon tracks, in upcoming neutrino telescopes, as the principal observable in 
our analysis. This ratio depends on the neutrino energy, density profile of the sources and on the 
oscillation parameters. The largely unknown flux normalization drops out of our calculation and 
only affects the statistics. For the current knowledge of the oscillation parameters we find that the 
matter-induced effects are non-negligible and the enhancement of the ratio from its vacuum value 
takes place in an energy range where the neutrino telescopes are the most sensitive. Quantifying 
the effect would be useful to learn about the astrophysics of the sources as well as the oscillation 
parameters. If the neutrino telescopes mostly detect diffuse neutrinos without identifying their 
sources, then any deviation of the measured flux ratios from the vacuum expectation values would 
be most naturally explained by a large population of hidden sources for which matter-induced 
neutrino oscillation effects are important. 

PACS numbers: 96.40.Tv, 14.60.Pq, 98.70.Rz, 98.70.Sa 



I. INTRODUCTION 

High-energy (> 100 GeV) neutrino emission has been 
predicted from several types of astrophysical sources 
such as active galactic nuclei (AGNs), gamma-ray bursts 
(GRBs), core-collapse supernovae (SNe), supernova rem- 
nants (SNRs), microquasars, etc. All these sources 
have been observed in electromagnetic wave bands rang- 
ing from radio to high-energy 7-rays. Low energy (~ 
10 MeV) thermal neutrinos have been detected from a 
nearby core-collapse SN 1987A [2j. High-energy neutri- 
nos may be detected from a much longer distance because 
of an increasing interaction cross-section with energy 
[USUI- Upcoming kilometer scale ice/water Cherenkov 
detectors such as IceCube @ in Antarctica and its coun- 
terpart in the Mediterranean, called KM3NeT @], will 
thus open up a new observation window in high-energy 
neutrinos. 

Neutrinos, unlike their electromagnetic counterparts, 
may carry information on temperature, density, etc. from 
deep inside the astrophysical sources. Their detection 
may also be used to probe neutrino flavor oscillation pa- 
rameters in matter and in vacuum as they propagate in- 
side the sources and over astrophysical distances to reach 
Earth. The results will be complementary to accelerator 
and reactor based experiments. Particular examples have 
been carried out in detail for ~ 10 MeV thermal neutri- 
nos from core-collapse SNe 0. 

Oscillations of high energy neutrinos have also been 
extensively discussed in Refs. d El El El ES EI 
EH EH E3- The sensitivity of high energy neutrinos 
to oscillation par ameters has been explored in Refs. 
0, EH El H US Matter effects in the context of 
neutrinos from dark matter annihilation have been con- 



sidered in Ref . pS] . Oscillations of solar atmosphere neu- 
trinos have been discussed in Ref. [lfj . Matter effects for 
high energy neutrinos are generally small. In order to 
get significant effects it is necessary to encounter a reso- 
nant density, as well as to g o through a minimum matter 
width, as shown in Ref. [17| • These conditions are usually 
not satisfied for optically thin sources. 

In this paper we explore flavor oscillation effects on 
high-energy neutrinos from optically thick sources where 
matter effects can be larger. These neutrinos are pro- 
duced deep inside the astrophysical sources (so-called 
"hidden sources") and as such their emission from the 
source is not accompanied by any high-energy electro- 
magnetic component. Examples of hidden sources may 
be highly relativistic GRB jets as they are formed inside 
collapsed stars and start burrowing their way through 
the stellar envelope [IS [lj| ; semi-relativistic jets inside 
core-collapse SNe which may impact and disrupt the en- 
velope but choke inside [20(; and core dominated AGNs 
21, 22]. The first two cases are motivated by the obser- 
vation of a number of GRBs associated with supernovae 
[47l | which supports the collapsar model of GRBs [23| . 
The collapsar model ties the GRBs and supernovae in a 
common thread, both originating from core-collapse of 
massive stars. While observed GRBs are endowed with 
highly relativistic jets, many more such collapses are ex- 
pected to produce mildly relativistic jets. In both cases 
the jet is launched from inside the star [IS [H, HH] and 
burrows through the stellar interior while the envelope is 
still intact. High energy neutrinos, produced by collisions 
of plasma materials (so called "internal shocks") in the 
jet, are emitted following a density gradient while the jet 
itself may or may not break through the envelope. High- 
energy neutrinos are produced via decays of pions and 
kaons created by hadronic (pp) and/or photo- hadronic 
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(pj) interactions of shock-accelerated protons in these 
examples. We concentrate on matter enhancement of the 
electron (anti)neutrino flux due to the small mixing an- 
gle #13 as these neutrinos propagate to the stellar surface 
from their production site. Matter effects have not been 
previously discussed for these sources and they provide 
an opportunity for identifying and studying such "hid- 
den" sources through the neutrinos they emit. We show 
in particular the energy range where this effect may be 
detectable by neutrino telescopes, and we also comment 
on possible extraction of neutrino properties such as mass 
hierarchy and the CP violation phase. 

The organization of the paper is as follows: We out- 
line the astrophysical source model and neutrino flux 
parametrization in Sec. [ill We discuss neutrino oscil- 
lations, our numerical approach to calculate the effects 
and analytic expectations in Sec. IHIl Detection and our 
results are given in Sec. IIVI and in Sec. [V] respectively. 
Conclusions and outlooks are given in Sec. EU 

II. ASTROPHYSICS OF THE SOURCE 

The particular hidden source model we employ for cal- 
culation purpose is a jetted core collapse supernova model 
in Ref. [2fj. Our treatment may however be applied to 
any generic case of a neutrino flux and matter density 
profile. The presupernova star is a blue supergiant (BSC) 
with a radiative hydrogen envelope which is capable to 
produce a Type II or lb SN. The hydrogen envelope sits 



on a helium core of radius trc ~ 10 11 cm. In case of 
significant stellar mass loss, the presupernova stars may 
have no hydrogen envelope left at all and produce Type 
Ic SNe, a fraction of which are now strongly believed to 
produce GRBs as well. However, we do not investigate 
this case as the relativistic jet front (rj e t ~ 10 10 8 cm) 
is too close to the stellar surface to have any significant 
amount of material in front [20| . 

An analytic form of the density distribution near the 
edge of each layer of a star with polytropic structure is 
given by p = pi(H k /r — 1)" [26]. Here R+ is the star's 
radius and the polytropic index n — 3 for a radiative 
envelope with constant Thomson opacity. For a helium 
star of _R* w 10 11 cm, pi w 2 g cm -3 and n = 3 were 
found by fitting data from SN 1998bw [28|. The presu- 
pernova star of SN 1987A is a BSG with R* » 3- 10 12 cm, 
pi f» 3- g cm~ 3 and n = 3 (H, [3(|. Below we write 
3 models of the density profile that we use for a BSG of 
i2* = 3 ■ 10 12 cm. All models are normalized to give the 
same density p = 2(10 n /10 10 - 8 — l) 3 = 0.4 g cm -3 , as in 
the case of a helium star, at r = rj ct = 10 10 8 cm. 

Model [A] corresponds to a polytropic hydrogen enve- 
lope with p(r) tx r -3 , scaling valid in the range rj ot < 
r < R+. Model [B] is a power-law fit with an effective 
polytropic index n c g = 17/7 as done for SN 1987A in 
Ref. [27|. Model [C] includes a sharp drop in density at 
the edge of the helium core. Note that our choice of dif- 
ferent n e ft below and above the helium core, in this case, 
is motivated by Ref. [26[ . The parameter A corresponds 
to the drop in density and its value is set by hand. 



[A] p(r) = 4.0 ■ 10~ 6 ( ^ - 1 ) & cm ~ 3 (!) 
™ , n o.^-5 f {RJr) 17/7 ; 10 10 - 8 cm < r < r b = 10 12 cm , 

[B] P(r) — 3.4- 10 X < M7/7/ „MS„ „Mi . § cm (2) 



[nl M en , n - 6A fR* 1 A"°' T -3 , (2.1,20) ; 10 10 - 8 cm< r < 10 11 cm 

[C] p{r) = 6.3 -10 M 1 g cm 3 ; (n cS ,A) = < ) n lnll (3) 



r ) v ' I (2.5,1) ; r > 10 11 cm 
I 

We have plotted the density profiles in Eqs. (TT]), and parametrized as [3J 

([3]) in Fig. [TJ The number density of electrons is given 

by N e (r) = N A p(r)Y e , where Y e m 1 in our cases, is the , / d L \~ z { E, 

$1 „ = 4 • 10 ~ 



number of electrons per nucleon or the electron fraction. T,"it ^ • ru \ cra J \ 300 GcV 

As for reference, N e — 2.4 • 10 23 cm~ 3 at r = rj ot = _ x _ 2 _ 1 

pQio.8 cm for all moc jels described above. GcV cm s 



100 < E v /GeV < 3 • 10 5 , (4) 



= 8 • 10 38 ' 



d r ,\ 2 ( E„ 



We use the pion and kaon decay neutrino flux mod- K ^ fi {cm J ^5 • 10 4 GeV / 

els from hadronic (jpp) interactions by shock acceler- GeV _1 cm _2 s _1 ■ 

ated protons in the semi-relativistic hidden jets [13, Hl[ . / 4 

The fluxes = <f>^ = 2$ Ve = 2$ S J, from an iso- K = J 3 5 100 ~ ^/ GcV ^5-10 

lated source, at the production site r = rj et may be 1 4 ; 5 • 10 4 < E^/GeV < 3 • 10 5 



(5) 
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FIG. 1: (Color online) Models of outer density profile of a 
blue supergiant star of radius R+ = 10 12 ' 5 cm. The outermost 
hydrogen envelope sits on a helium envelope extending up to 
a radius of ~ 10 11 cm. High-energy neutrinos are produced at 
a radius rjet = 10 10 ' 8 cm by a relativistic jet. The models [A], 
[B] and [C] are described in Eqs. fT}, @ and © respectively. 



respectively from pion and kaon decays. Here cLl is the 
luminosity distance of the source. Next we discuss oscil- 
lation effects on these neutrinos as they propagate from 
inside the stellar interior to the surface in dense me- 
dia, from the source to Earth and through the Earth. 
Note that high energy neutrinos are produced in shocked 
material in the jet which has a density lower than the 
surrounding stellar material. However, the width of the 
shock is too small (jet radius divided by the Lorentz boost 
factor of the jet) to have any significant oscillation effect 
before the neutrinos start moving through the jet head 
(which is roughly in equilibrium with surrounding mate- 
rial) and the envelope. Also, the jet is well outside the 
stellar core which may be turbulent and the envelope is 
not disrupted yet. 



III. NEUTRINO OSCILLATION IN VACUUM 
AND IN MATTER 

Solar, atmospheric, accelerator and reactor neutrinos 
have provided ample evidence for neutrino oscillations. 
Neutrinos from astrophysical sources like those described 
above are affected by oscillations, which change the fla- 
vor composition of the fluxes between production and 
detection. 

The neutrino flavor eigenstates u a , where a = e, /i, r in 
the case of 3 flavor mixing, are related to the mass eigen- 



states Vj, where j = 1,2,3 corresponding to the masses 
in j , by the Maki-Nakagawa-Sakata (MNS) unitary mix- 
ing matrix U as v a — U a jVj. The sum over repeated 
indices is implied. We use the standard expression of U 
from Ref. [33| . 

We use standard oscillation parameters obtained from 
global fits to neutrino oscillations data: 

Ato 2 o1 = Amar - 8 x 1CT 5 eV 2 
#soi = 6 12 ~ 33.83° 
\Aml tm \ = |Am 2 2 |~2.4xKr 3 eV 2 

#atm = #23 ~ 45° 

sin 2 (26»i 3 ) < 0.15 (6) 

The sign of the atmospheric mass difference has not been 
determined. Positive (negative) Am 2 tm correspond to 
the normal (inverted) hierarchies. 

Uncertainties in the oscillation parameters are still rel- 
atively large and small variations due to these uncertain- 
ties on the generic results we present here are to be ex- 
pected. We will explore the impact of the expected er- 
rors of the oscillation parameters on our observable and 
we will discuss how these errors could affect an eventual 
measurement of the neutrino mass hierarchy (i.e normal 
vs inverted) via matter effects. 

We will also consider the effects of a CP violating phase 
which is at present unconstrained. A precise approxima- 
tion (for constant density and energies high enough such 
that A 2 i/V e <C 1 and A 2 iL Cl) for the v e — ► Vy, transi- 
tion probability is given by [34] : 

sin 2 6 23 sin 2 26 13 (^) * sin 2 (%) 

+ cos 2 23 sin 2 20 12 (Aai) 2 s in 2 (Y) (7) 
+J^^sin(Y) S in(¥)cos(±<5-^) , 

where Ay = Amf j /2E l/ , V T = \V e T &3i\, the matter 
potential V e = V e (r), J is the Jarslog invariant and the 
sign minus (plus) refers to neutrinos (antineutrinos). 

For the above values of Am 2 , even for the high energies 
of interest here, the distances relevant for the sources de- 
scribed above are so large that, after propagation through 
vacuum, only the averaged oscillation is observable for 
the simple two flavor case with the probability given by 

P{u a - „„) = sin 2 (2tf) sin 2 {^-) ^ \ sin 2 (20) . 

(8) 

In the three flavor case, given fluxes at the surface of the 
source cj/^ , <jf>°, and cpP, T , the fluxes at Earth are given by: 

K =€*-\ sin 2 2012(20°. " <f>% ~ O (9) 
K =^r=\^% + 4>l r ) + \ sin 2 20 12 (2< - 4>l. - O 
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if we assume that #23 is maximal and #13 is very small (as 
indicated by neutrino oscillation data). If neutrinos are 
produced in pion and/or kaon decays, the initial flavor 
ratio is given by v e : : v T = 1 : 2 : 0. After propagation 
over very long distances in vacuum, neutrino oscillations 
change this ratio to 1 : 1 : 1 because of the maximal 
Vp v T mixing. For the sources we are considering, in 
the low energy range, even the distance traveled inside 
the source is large enough that neutrinos oscillate many 
times and the phase information is lost. When the initial 
ratio v e : z/„ : v T is different from 1:2:0, the flavor ratio 
at Earth is affected by the full three flavor mixing and is 
different from 1:1:1. 

In the case of interest here, even though neutrinos are 
produced by pion and kaon decays, the fluxes at Earth 
can have a flavor composition quite different from the 
standard 1:1:1 ratio. This is because the neutrino 
flavor composition is affected by the propagation through 
matter inside the astrophysical object. 



The fluxes at the surface of the star are given by the 
sum of products of the fluxes at the production site 
[Eqs.Q and ©] and the oscillation probabilities as 
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(10) 



The same type of relations apply to anti-neutrinos. 

The v a <-> vp oscillation probabilities P a( g on the stellar 
surface r = R+ can be obtained by solving, with the 
appropriate initial conditions, the evolution equation: 
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Here V e (r) — \/2GpNe{r) is the matter-induced charged 
current potential for v e in a medium of electron number 
density N e (r) found from Eqs. © and ©. For v e 
the potential changes sign. 

Matter effects are expected to be large when neutrinos 
go through a resonant density. For constant electron den- 
sity, there would be two such resonances, corresponding 
to the two mass squared differences: 



N? 



Am 2 ol cos 26*12 
2V2G F E 



I Am? 



I cos 29 
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2V2G F E 



3 x 10 23 

E\GeV] CUi 
^ 9 x 10 24 
~ E[GeV] ' 



(12) 
(13) 



As can be seen from the previous section, these densities 
can be encountered in the sources of interest by neutrinos 
with the relevant energies. 

Even if #13 is very small, matter effects can be expected 
from the solar Am 2 oscillations [see Eq. (p~2|) ] . If 6*13 is 
non-negligible, larger matter effects can be expected from 
the atmospheric Ato 2 oscillations [see Eq. (fT3"|) ]. These 
are non-negligible at a relatively higher energy. 

Since the matter potential is different for anti- 
neutrinos, it is important to independently solve the evo- 
lution equations for neutrinos and anti-neutrinos and 
consider the averaging over the two fluxes only when 
reaching the detector, which cannot distinguish between 
the two. Note that if the hierarchy is normal (inverted) 
the resonant behavior occurs only in the neutrino (anti- 
neutrino) channel. 



The matter density in the sources we are considering 
is varying rather strongly, such that a resonance is not 
directly observable. Adiabatic conditions are not usu- 
ally satisfied, such that a full numerical treatment of the 
problem is necessary. The adiabatic approximation can 
still be used in some energy range. It is useful to under- 
stand some of the features introduced by matter effects 
by analyzing limiting regions of the models described in 
the previous section. The adiabaticity parameter is de- 
fined as (see Ref. |8|] and references therein): 



Am 2 sin 2 29 



1 



7 



2E cos26> (l/N e )(dN e /dr) 



(14) 



i.e. the ratio of the resonance width and the neutrino 
oscillation length. In other words, the adiabaticity pa- 
rameter represents the number of oscillations that occur 
in the resonance region. Adiabatic transition (7 3> 1) 
means that the neutrino will remain in its particular su- 
perposition of initial instantaneous mass eigenstates as 
it crosses the resonance. The "flip" probability, that is, 
the probability that a neutrino in one matter eigenstate 
jumps to the other matter eigenstate is 



( n 

Ffiip = exp \-- 



(15) 



which means that, if the neutrino adiabaticity param- 
eter is 7 ^> 1, the "flip" probability Pfl ip ~ 0, many 
oscillations will take place at the resonance region and 
strong flavor transition will occur. On the other hand, 
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if the adiabaticity parameter is 7 > 1, only a few oscil- 
lations take place in the resonant region and the mixing 
between instantaneous matter eigenstates is important, 
and Pflip gives the jumping probability from one instan- 
taneous mass eigenstate to another. For the solar tran- 
sition, where the mixing angle is large, the appropriate 
expressions for the adiabaticity conditions and flip prob- 
abilities can be found in [35j ■ Since the main effects that 
we are interested in correspond to the small mixing angle, 
we limit our analitical discussion to this case. We would 
like to emphasize that in order to obtain the correct neu- 
trino oscillation probabilities a full numerical treatment 
of the problem is necessary and our discussion based on 
flip probability serves only as a qualitative guide to un- 
derstanding the solution in some limited energy ranges 
where adiabaticity is satisfied. 

Figure [2] shows the flip probability corresponding to 
the atmospheric transition for the three density profiles 
we consider. We explicitly delimit three different regions: 
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FIG. 2: (Color online) Flip probability corresponding to the 
atmospheric transition for sin 2 2613 — 0.15, for the three den- 
sity profiles in Fig. [T] 

the adiabatic region corresponds to very small flip prob- 
abilities (Pflip < 0.1) and strong flavor transitions; in the 
intermediate region (0.1 < Pfli P < 0.9) the transitions 
are not complete; the highly non-adiabatic region corre- 
sponds to very high flip probabilities (Pflip > 0.9). 

At low energies transitions are adiabatic for the density 
profiles in models [A] and [B] , so strong conversion is ex- 
pected for these two models. Model [C] has a sharp drop 
in density and at energies below 500 GeV the resonant 
density is reached only on the step, where it is highly 
non-adiabatic, so one does not expect a significant mat- 
ter enhancement. 



At higher energies the adiabatic regime applies only 
for Model [B], for which oscillation effects are expected 
to be large, while for models [A] and [C] the flavor transi- 
tion is incomplete and a smaller effect is expected. These 
expectations are confirmed by our exact results, that is, 
around energies ~ 1 TeV, the effect for the Model [A] 
and the Model [C] should be similar and smaller than 
the one observed for the Model [B]. The differences in 
the adiabaticity behavior for the three models are in- 
duced by their different matter profiles, see Fig. [TJ for 
Model [B] (which has the least steep matter density), the 
adiabaticity condition is satisfied for a larger neutrino en- 
ergy range. The adiabaticity of the transitions depends 
on the matter density profile. If this is very steep, the 
corresponding resonance width will be very small (it is 
inversely proportional to the derivative of the matter po- 
tential) and consequently the neutrino system can not 
adapt itself, only a few oscillations will take place and 
the flavor transitions will not be complete. In order to 
have large matter effects it is also necessary to go through 
a minimum matter width, as shown in Ref. [17j | . For the 
sources we consider here the density is high enough that 
this condition is satisfied. 

In Fig. [3] we show the electron, muon and tau neutrino 
fluxes at the surface of the source, normalized to the ini- 
tial electron neutrino fluxes. We compare the case when 
only the vacuum oscillations are considered inside the 
source with the fluxes of neutrinos when matter effects 
are taken into account with a density profile as in Model 
[A] or Model [B], for normal hierarchy. The averaging 
due to fast vacuum oscillations can already be observed 
in the lower energy range. It can also clearly be seen that 
matter effects modify the flavor composition of the neu- 
trino fluxes, introducing energy-dependent features. As 
expected from the discussion above, the effects are larger 
for Model [B] where adiabatic-like transitions occur in 
a wider energy range. We also show the corresponding 
anti-neutrino fluxes for Model [A] for comparison. The 
results are very close to those in vacuum, as expected for 
the normal hierarchy. 

After propagating the neutrinos through matter to the 
surface of the star, the fluxes at Earth immediately fol- 
low by considering the averaged oscillations over the long 
travel distance in vacuum as obtained from Eq. (|10p . 

Figure 2] makes the same comparison presented in Fig. 
[3] after propagation all the way to the Earth. The flavor 
composition of the neutrinos gets further modified by the 
propagation from the source to the Earth. 



IV. DETECTION 

On their way to detectors, neutrinos also propagate 
through matter in the Earth. In this case matter effects 
on neutrino oscillations are however extremely small, 
since the energies considered here are much higher than 
the resonant energy inside the Earth. For Am^ tm = 
2.4 x 10 -3 eV 2 the corresponding resonant energy for 
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FIG. 3: (Color online) Neutrino fluxes at the surface of the 
source (for sin 2 2#i3 = 0.15), normalized to the initial electron 
neutrino fluxes. From upper to lower panels: vacuum, Model 
[A], Model [B] and Model [A] antineutrinos. The solid, dashed 
and dotted curves are for v e , and v T respectively. 



a neutrino going through the Earth's mantle (with the 
density ~ 3 g/cm 3 ) is around 10 GeV, already lower 
than the energies we consider. The resonant energy is 
even lower for the solar mass difference or for trajecto- 
ries going through the higher density Earth core [3q |. 

The charged current and neutral current neutrino- 
nucleon interaction cross-section increases with energy 
and attenuation effects start becoming important for the 
propagation through matter. The interaction length: 



C 



1 



o vN (E)N A 



(16) 



becomes comparable to the Earth diameter around 40 
TeV. The attenuation effects are thus becoming relevant 
only at the highest energies considered here, where the 
fluxes are very small. Consequently, propagation through 
the Earth will have negligible effects on the neutrino 
fluxes and flavor composition previously discussed. 

Neutrinos from the sources discussed here are expected 
to be detectable in IceCube and other neutrino telescopes 
which have good sensitivity in the energy range of inter- 
est, between 100 GeV and 100 TeV. 

It is possible in IceCube to get information about the 
flavor composition of the neutrinos as well. The detector 
is mostly sensitive to observing muons identified by their 
very long tracks, thus counting the number of interac- 
tions. Electron neutrino interactions create electromag- 
netic cascades which can be observed in IceCube. Tau 
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FIG. 4: (Color online) Same as Fig. [3] but at the surface of 
the Earth. 



neutrinos will also lead to cascades. At these energies the 
tau decay length is very small and the interaction of the 
neutrino and tau decay cannot be separated. It might be 
possible to separately infer the fluxes of all three flavors 
if hadronic showers can be well separated from electro- 
magnetic showers through their muon track content. We 
will consider taus to be indistinguishable from electrons 
in our analysis and we will compare the number of tracks 
due to and number of showers due to v e and v T . 

Another important feature of interest for the detection 
of the effects studied in this paper is the energy depen- 
dence of the signal: as discussed in the previous section, 
matter effects have a strong dependence on energy, which 
is correlated with the density profile inside the source; 
changes in neutrino oscillation parameters also induce 
energy dependent effects, as we discuss in more detail in 
the next section. The energy resolution in IceCube is ex- 
pected to be somewhat better than 30% on a logarithmic 
scale for muon neutrinos and about 20% on a linear scale 
for cascades. 

We explore the spectral shape of the shower-to-muon 
track ratio inferred at Icecube, defined as: 



R 



(17) 



The absolute flux normalization thus does not affect our 
results. The number of events, N Vot i 0a \ is proportional 
to: 



l a (p a ) +< <p a ) ) (E v ) x u Va (Pa )(E v )dE v , 

(18) 
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i.e. we have computed the (anti)neutrino fluxes from pion 
and kaon decays at the detector after considering the 
propagation inside the source (affected by matter effects) 
as well as the propagation from the source to the detec- 
tor. The (anti)neutrino fluxes are then convoluted with 
the (anti)neutrino cross-sections [4| and integrated over 
the neutrino energy, assuming a conservative energy bin 
size of AE V = 0.3E^. If the matter potential inside the 
astrophysical source is neglected, the neutrino flavor ra- 
tio in Eq. (jTTJ) should be almost constant, R ~ 2, over 
all the neutrino energy range. Adding the matter contri- 
bution to the neutrino propagation will change this con- 
stant value of R in an energy dependent fashion. In the 
following section we explore the spectral behavior of the 
shower-to-track ratio in the 300 GeV- 300 TeV neutrino 
energy range, considering both normal and inverted neu- 
trino mass orderings and exploring the effects of a non 
zero leptonic CP violating phase 8 and of uncertainties 
in other oscillation parameters. 



V. SPECTRAL ANALYSIS 

We present here the main results of our study. We are 
mainly interested in neutrino propagation through mat- 
ter inside the source. This is followed by the propagation 
through vacuum between the source and Earth. Matter 
effects and absorption inside the Earth are negligible for 
the energies considered here. As discussed in the previous 
section, a good observable for studying the effects of neu- 
trino propagation inside the source is the spectral shape 
defined in Eq. (17]). While R ~ 2 for the optically thin 
sources, energy dependent deviations from this value are 
expected due to matter effects in the astrophysical hid- 
den sources considered here. 

Figure [5] shows the shower-to-muon track ratio in 
Eq. (fT7|) for the normal and inverted neutrino mass hi- 
erarchies for the Model [A] density profile [see Eq. ([Tjl], 
Deviations from the standard scenario are non-negligible 
in the energy range around a few TeV. The impact of 
the matter potential is much larger if nature has chosen 
the normal mass ordering, since, if that is the case, only 
neutrinos can go through the resonance density and 
the effect will be larger due to the higher neutrino cross- 
section, when compared to the antineutrino one. For 
neutrino energies E v < 6 TeV, where the resonant effect 
in the (anti)neutrino propagation is expected to be lo- 
cated, the neutrino cross-section is roughly twice the an- 
tineutrino cross section. Therefore, the matter potential 
impact in the inverted hierarchy situation (when only the 
antineutrino propagation gets distorted) is much smaller 
than in the normal hierarchy case. We have also studied 
the shower-to-muon track ratio spectral shape when the 
CP violating phase S is not zero. The different curves in 
Fig. [5] correspond to different values of the CP violating 
phase S. They were obtained by numerically computing 
the full oscillation probabilities, but the relative shape of 
these curves with respect to the 5 = curve shape can be 
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FIG. 5: (Color online) The upper (lower) panel depicts the 
shower-to-muon track ratio defined in Eq. (|17|) for sin 2 2#i3 = 
0.15, for Model [A], for the case of normal (inverted) hierarchy 
and for different values of the CP violating phase 8. 



easily explained by means of the oscillation probabilities 
in Eq. ([7]). This is because the expression applies for 
almost all the relevant energies and baselines explored in 
our study. It can be easily seen, for example, that for 
a fixed distance of L = 10 5 km and a (constant) density 
correspnding to the electron number density for Model A, 
Ai 2 /V e < 1 and A 12 L < 1 for energies ~ 10 3 GeV and 
the approximated probability formula provides an accu- 
rate description of the oscillation probabilities versus the 
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CP phase 8. 

There are three very different terms in Eq. j7]): the 
first one, which is the dominant one, is responsible for the 
matter effects; the second one, which is the solar term, is 
dominant for very small values of #13 not explored here, 
and therefore, negligible for the present discussion; the 
third term, (named the interference term in the litera- 
ture), is the only one which depends on the CP violat- 
ing phase (5 and it determines the shape of the curves in 
Fig. In the case of neutrino transitions and normal hi- 
erarchy, notice that if 8 = 0, the three oscillatory factors 
of the third term can be written as sin (V1L/2). If, for 
instance, 8 — w, the three oscillatory functions can be 
written as — sin (V-L/2), which explains why when the 
8 = curve has a maximum and the 8 = tt curve has a 
minimum. If 8 = 7r/2, the three oscillatory factors can 
be simplified as sin (VLL/2) cos(V r _L/2), which explains 
the relative shift in the phase of the 8 — ir/2 curve with 
respect to the 8 = 0, it curves. Finally, if 8 = 3ir/2, the 
oscillatory dependence goes as — sin(V_L/2) cos(VLi/2), 
i.e, it is the opposite to 8 = ir/2. In summary, one 
can predict the shape dependence on the CP phase of 
the shower-to-muon track ratio just by observing the dif- 
ference in the oscillation pattern of the four functions 
± sin 2 x, ± sin(x) cos(ic). 

In the case of the inverted hierarchy, see Fig. \E\ (lower 
panel), the situation is reversed: notice, [see Eq. (JJJ)], 
that a change in the sign of Am^ can be traded in the 
vacuum limit (V e — > 0) for the substitution 8 ^ -k — 8. 
Matter effects, however, can break this degeneracy. By 
making use of the matter effects, in principle, it could be 
possible to extract the neutrino mass hierarchy, since the 
neutrino oscillation probability is enhanced (depleted) for 
positive (negative) Amf^. 

Figures [5] and [7] depict the shower-to-muon track ratio 
for the density profiles given by the models [B] and [C] 
in Eqs. ([2]) and ([3]) respectively; and the previous discu- 
ussions on Model [A] in Fig. \5\ apply to these cases as 
well. 

A change in the dependence of the matter potential 
profile will affect the location of the energy at which the 
maximum matter effect is located. Depending on the 
specific form of the electron number density vs distance 
N e (r), a neutrino with a fixed energy E v will reach the 
resonant density at a different distance from the center 
of the astrophysical source, as can be seen from Fig. [5J 
where we have depicted the distance at which the res- 
onant density is reached vs the neutrino energy E v for 
the three models. Notice that, for Model [C], there is a 
range of energies (between 100 and 500 GeV) for which 
the neutrino can profile, where the transition is highly 
non-adiabatic. The propagation inside the source from 
the production point to its surface is thus different mea- 
sured spectral shape. 

It would be interesting to explore if it is possible to use 
neutrinos from astrophysical sources as those discussed 
here in order to explore neutrino properties like the neu- 
trino mass ordering or the CP violating phase 5. 
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FIG. 6: (Color online) The upper (lower) panel depicts the 
shower-to-muon track ratio for sin 2 2#i3 = 0.15, for Model 
[B], for the case of normal (inverted) hierarchy, for different 
values of the CP violating phase S. 



A deviation from the expected 1:1:1 ratio has 
been extensively explored in the literature (sometimes 
involving highly exotic scenarios), e.g., a different flavor 
ratio at the source [H, [H, 53], decaying neutrino sce- 
narios [4lj],Pseudo-Dirac schemes [HI, additional sterile 
neutrinos (3| , magnetic moment transitions [45[ and the 
possibility of measuring deviations from maximal atmo- 
spheric mixing '12J. An extensive discussion of possible 
deviations from the equal flavours scenario, including os- 
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FIG. 7: (Color online) The upper (lower) panel depicts the 
shower-to-muon track ratio for sin 2 2#i3 = 0.15, for Model 
[C], for the case of normal (inverted) hierarchy, for different 
values of the CP violating phase 5. 



cillations and decays beyond the idealized mixing case, 
characterization of the source and CPT violation was pre- 
sented in [13]. Some of these studies predict energy de- 
pendent ratios, however, none discussed matter effects 
on high-energy neutrino flux ratios as we explore in this 
paper. 

If a relatively large deviation (~ 10%) of R from its 
standard value R ~ 2 is observed for energies E < 2 
TeV, one could infer a non zero value for the mixing angle 



FIG. 8: (Color online) Distance from the center of the as- 
trophysical source (in cm), at which the resonant density is 
crossed vs the neutrino energy, for Model [A] (solid line), 
Model [B] (dashed line) and Model [C] (long-dashed line), 
for sin 2 20 13 = 0.15. 



6*13. The energy location of the matter enhanced peaks 
in the shower-to-muon track ratio depends on the precise 
value of 6*13: for smaller values of sin 2 26>i3 < 0.15, the 
effect is larger (smaller) at lower (higher) energies with 
respect to the effect illustrated here for sin 2 26>i3 = 0.15, 
as expected from Eq. (fl~3| . Figure [9] shows the change 
in R when varying #13, for the Model A density profile, 
normal hierarchy and 8 = 0. It can be seen that the 
results strongly depend on the #13 mixing angle, whose 
effects get enhanced by a matter resonance for values 
larger than about sin 2 26>i3 = 0.01, when the transitions 
are mostly adiabatic. 

Measuring the CP phase by observing the shower-to- 
muon track ratio is highly challenging, since for that pur- 
pose a precise knowledge of the matter density profile 
p(r) would be required (see Ref. |46l| fo r the prospects 
in a decaying neutrino scenario, Ref. [10[ for a combined 
analysis of neutrinos from reactors and optically thin as- 
trophysical sources, and Ref. 

@, ElIIl El for astrophys- 
ical neutrinos in different scenarios). 

Under the assumption of a quite good knowledge of the 
density profile, it could be in principle possible to deter- 
mine/verify the neutrino mass hierarchy since in all three 
models the effect is much higher in the normal mass or- 
dering than in the inverted one regardless the value of the 
CP violating phase S. However, in practice, the errors on 
the remaining oscillation parameters might compromise 
the possibility of extracting the neutrino mass hierarchy, 
in particular the uncertainties on the atmospheric mixing 
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FIG. 9: (Color online) Variation of R with 9i 



FIG. 10: (Color online) Variation of R with #23. 



angle #23 (if a more precise knowledge of this mixing angle 
is not available at the time of the shower-to-muon track 
ratio measurement) . FigureHOlshows the effects of chang- 
ing 623 between maximal mixing and sin 2 #23 = 0.38 (the 
current 2 sigma allowed region) , for the Model [A] density 
profile, normal hierarchy and 8 = 0. Matter effects are, 
in general, larger for larger values of sin 2 623 , except for 
very high energies when the ratios approach their values 
in vacuum. In this range, if 823 is larger than 7r/4 the 
muon flux is enhanced and R < 2, while the opposite is 
true if 623 is smaller than 7r/4. The uncertainty in the 
mixing parameters would affect the possible determina- 
tion of the mass hierarchy since the inverted hierarchy 
pattern could be confused with a smaller value of 823. 
A separate measurement of electron and tau neutrinos 
would greatly improve this situation, allowing to extract 
both the hierarchy and the value of #23- 



VI. CONCLUSIONS AND OUTLOOK 

We have discussed here high-energy neutrinos pro- 
duced in optically thick (hidden) astrophysical objects for 
which large matter density inside the source can affect the 
oscillations of these neutrinos. The matter-induced tran- 
sitions can lead to significant deviations from the 1:1:1 
flavor ratios expected in standard scenarios. These devi- 
ations are expected at specific energies determined by the 
density profiles inside the sources. IceCube would be in 
an ideal position to measure such effects. The main ob- 
servable that would be sensitive to these effects and that 
we have studied in detail is the shower-to-muon track 
ratio defined in Eq. (fTT|) . 



A large number of events would be necessary in order 
to measure the effects of neutrino oscillations inside the 
source. In order to establish a 3 sigma effect, more than 
1000 events would be needed. This would require a rel- 
atively nearby source, within a few megaparsecs. Given 
the supernova rate in nearby starburst galaxies such as 
M82 (3.2 Mpc) and NGC253 (2.5 Mpc) is about 0.1/yr, 
the possibility of such an event to take place is not rare 
assu ming a significant fraction of SNe are endowed with 
jets [20l I3TL l32|. Upcoming neutrino telescopes will be 
able to constrain this fraction. The combined SN rate 
from all galaxies within 20 Mpc is more than 1/yr [37j ■ 

With such large number of events it would also be pos- 
sible, in principle, to investigate neutrino properties like 
the mass ordering and CP violating phase, especially if 
future reactor and accelerator experiments could provide 
a better measurement of the mixing angles. If all neu- 
trino properties would already be known from other ex- 
periments, the shower-to-muon track ratio measured in 
IceCube could be used to determine source properties: 
matter effects inside the source reflected in the ratio R 
are extremely sensitive to the source density profile. 

Since high-energy 7-rays are not emitted from opti- 
cally thick sources, it is hard to predict their occurrence 
rate. If a large population of hidden high-energy neu- 
trino sources would exist in nature, then their combined 
effect would be evident in the ratio R measured from dif- 
fuse fluxes. A modulation of R with energy over a wide 
range is expected if the astrophysics varies from source 
to source and/or if the sources are distributed over a 
wide redshift. This is a distinct signal from the effect, 
e.g., discussed in Refs. (38l . |39( |: which appear only at 
ultrahigh-energies. 
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Neutrino telescopes are sensitive to the atmospheric 
neutrinos as well in the 0.1-100 TeV energy range we are 
interested in. However, the atmospheric neutrinos should 
not affect measurement of the ratio R for a nearby tran- 
sient point source that we are considering. For diffuse 
flux from all such point sources, the measured ratio will 
be a convolution of the atmospheric flux ratios and the 
astrophysical flux ratios. The oscillation effects for the 
atmospheric neutrinos is rather small (< 10%) because 
of a small Am 2 L/E. The experiments should be able to 
de-convolute them even without precise knowledge of the 
oscillation parameters. Assuming neutrino oscillation pa- 
rameters are well-measured, matter oscillation effects as 



we discussed here would be the most natural explanation 
of any deviation of R from its value in vacuum measured 
by upcoming neutrino telescopes. 
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